function PIE=realPIE(PIE,S,i);
  ss=S(i,:);
  x=find(sum(PIE>ss,2)>0);
  PIE(x,:)=[];
end
